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Finite temperature properties of symmetric ±J multileg Ising ladders and tubes are investigated using 
the statistical transfer matrix method. The temperature dependences of the specific heat and entropy are 
calculated. In the case of tubes, it is found that the ground-state entropy shows an even-odd oscillation 
with respect to the number of legs. The same type of oscillation is also found in the ground-state energy. 
On the contrary, these oscillations do not take place in ladders. From the temperature dependence of the 
specific heat, it is found that the lowest excitation energy is 4J for even-leg ladders while it is 2J otherwise. 
The physical origin of these behaviors is discussed based on the structure of excitations. 
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1. Introduction 

The spin systems with ladder and tube geometries 
have been attracting the interest from various view- 
points. 1,2 In the case of the spin-1/2 quantum Heisenberg 
ladders, it is well known that the properties of the ground 
states of the even-leg and odd-leg ladders are essentially 
different. Recently, the spin tube materials are also syn- 
thesized 3 and activated their theoretical studies. Various 
exotic quantum phenomena arising from the interplay of 
quantum fluctuation and frustration are predicted. 

Although the quantum spin ladders and tubes are ex- 
tensively studied, their classical counterparts have been 
less studied. Actually, the ground states of the regular 
unfrustrated classical Ising ladders and tubes are rather 
trivial. However, the ground state in the presence of the 
quenched randomness and frustration is nontrivial even 
in the Ising models. 4 ~ 6 Among them, the two-leg ±J 
Ising ladder is one of the simplest models with random- 
ness and frustration. Mattis and Paul 4 (MP) proposed 
the method to calculate the free energy of this model 
exactly, using the statistical transfer matrix method. Al- 
though their contribution was pioneering, their argument 
was limited to the case of two-leg ladder, and the numeri- 
cal estimation of the free energy and ground-state energy 
was inaccurate. In this paper, we extend their method to 
the multileg ladders and tubes to calculate their free en- 
ergy, entropy and specific heat at finite temperatures. At 
low temperatures, it is shown that the entropy of the 
classical ±J Ising tubes shows an even-odd oscillation 
with respect to the number of legs q, while the entropy 
of the classical ± J Ising ladder shows no oscillation. It is 
also found that the energy gap is 4 J for even- leg ladders 
but it is 2 J in tubes and odd-leg ladders. These features 
are qualitatively understood by considering the structure 
of the excitations. 

The even-odd oscillation in the correlation length of 
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the ±J Ising tubes was mentioned in ref. 7. However, 
these authors were interested in the limit of two dimen- 
sional ±J model (q — > oo) and this effect was consid- 
ered as a finite size effect which is harmful in taking the 
thermodynamic limit. Considering the recent increase of 
the interest in the models and materials with ladder and 
tube geometries, however, these peculiar properties of 
±J Ising ladders and tubes should be investigated in 
more detail. 

The present paper is organized as follows. In the next 
section, we introduce the model Hamiltonian. The sta- 
tistical transfer matrix approach by MP is extended to 
the multileg ladders and tubes in §3. The numerical re- 
sults are presented in §4. The last section is devoted to 
summary and discussion. We also present the correct es- 
timation of the ground-state energy and free energy of 
the two-leg ladder in Appendix. 

2. Models 

We consider the symmetric ± J Ising g-leg ladder 



H 



ladder 
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where S n . a (= ±1) is the Ising spin variable on the n-th 
rung and ev-th leg. The number of the rungs is denoted by 
L. The exchange constants J± na and J na are quenched 
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random variables which take the values ±J(J > 0) with 
equal probability. 

3. Statistical Transfer Matrix Formulation 

Before constructing the transfer matrices, we gauge 
out the randomness along the legs from the Hamiltonians 
(1) and (2) by the transformation 



respectively, where t n , a 's(= ±1) are quenched random 
variables. 

In this representation, the trace over T n in the par- 
tition function Z({t na }) can be readily taken for each 
realization of {t n>a }. This yields 



Z({t n , a })=Tr L l[V({t n , a }), 



(12) 



"no — I II s gnJ n ' Q J b n 



(3) 



W=l 



The Hamiltonians (1) and (2) are transformed into the 
forms, 

rrladdcr \ "* \ "* j c Q 

n— 1 a — 1 



where V({t na }) are 2 9_1 x 2 q ~ 1 sized transfer matrices 
between neighboring rungs parameterized by {t na }. In 
the following, we index the state of a rung {L a } by i = 

1 + Ea=i((^« + !)/ 2 x 2 " _1 )- Then > the elements of V 
are given by 



^.r cr ({U) = exp 



T 



L q 
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respectively, with 

J±na = J±na J\ SgnJ„< Q JT SgnJ„', Q+ i = ±J. 



\n' = l 



\n' = l 



(6) 



We define the spin variables T n and L na by 

T n = Sn.iSn+i^i, (7) 

Lna = S n . a S n ,a+l, ((X = 1, ...,q — 1). (8) 

It also follows that 

9-1 _ 

S n ,lS n ,q =J_J_ L na . (9) 

a=l 

Using these relations, the Hamiltonians (4) and (5) are 
further transformed into the forms, 

L <j-l 

rrladdcr j \ " \ " . f 
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for ladders, and 

v^ bc ({U) = exp 1 3 - r£ t a i a + *, j] z Q j I 

x2co S h|^h+^f[^L'j| (14) 

for tubes, where T is the temperature. 

Denoting the statistical mechanical weight of the i- 
th state of the n-th rung by x n j, we define the weight 
vector on the n-th rung x n = {x n _\, x n ,2, •■, XnZQ- 1 ) with 

2 q-l 

normalization > x„,j = 1. Then, the weight vector on 

i=l 

the (n + l)-th rung is determined by 

9 q-l 
^ 1 

%n+l,j = V^ X n ,iVij ({t a }) (15) 

Xn+i fri 

where Xn+i is the normalization constant for x n+ i de- 
termined by 

2 q-l 2 q-l 

Xn+l(x n , {tn,a}) = J^ ar n ,i J^ V it j ({t n ,a})- (16) 

i=l j=l 

It should be noted that 2« _1 (2 q ) vectors x n+x and 2 9 " 1 
(2 9 ) scalars Xn+i are generated from a single vector x n 
by the 29- 1 (2« ) possible choices of {i„, Q } in V" laddcr 

(ytube y 

Following MP, the free energy per spin F/N is ex- 
pressed using x„'s as 



(17) 



n=l 



where (...) , t -, means the average over {t n ^ a } and 7V(= 
Lq) is the number of spins. 
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Fig. 1. Temperature dependence of the specific heat per spin for 
(a) ladders and (b) tubes. The inset of (b) shows the temperature 
dependence of the latter in the low temperature regime. 

(a) 




Fig. 2. Temperature dependence of the entropy per spin for (a) 
ladders and (b) tubes. 




Fig. 3. Ground-state entropy of ladders (•) and tubes (o) per 
spin plotted against q. 




Fig. 4. Ground-state energy of ladders (•) and tubes (o) per spin 
plotted against q. 




Fig. 5. Ground-state entropy of ladders (•) and tubes (o) per 
spin plotted against 1/q. The right-directed open triangle is the 
value for the two-dimensional ± J Ising model. 
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Fig. 6. Ground-state energy of ladders (•) and tubes (o) per spin 
plotted against 1/q. The right-directed open triangle is the value 
for the two-dimensional ± J Ising model. 



4. Numerical Results 

For the numerical calculation, we fix the weight of the 
boundary state as x\_ a — l/2 9_1 (a = 1,...,2 9_1 ) and 
generate x n and Xn iterating (15) and (16). After n — 1 
iterations we have 2( 9-1 )(™ -1 ) (2 9 (™~ 1 )) weight vectors x n 
on the n-th rung for ladders (tubes). For large enough 
n, they correspond to the bulk weight and the effect of 
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Fig. 7. Arrhcnius plot of the low temperature specific heat of (a) 
ladders and (b) tubes. In the inset of (a), the interval 2.5 < J/T < 3 
is magnified . 



the fixed boundary weight is washed out. Taking into ac- 
count the self-averaging nature of ln^n, we obtain the 
thermodynamic limit of F/N by averaging — (T/q)lnx n 
over all possible Z^- 1 )^- 1 ) (2«(™- 1 )) values of %„. With 
the increase of q, however, the summation over all pos- 
sible \t n ^a\ becomes too demanding. Therefore, in our 
calculation, the average is taken over randomly chosen 
1000 realizations of {t n . a } with 200 < n < 1200. 

The specific heat and entropy are calculated by the 
numerical differentiation of the free energy. We calculate 
the free energy at different temperatures for the same set 
of {t n , a }- Since the free energy for each set of {t n , a } is 
a smooth function of temperature, the averaged free en- 
ergy is also a smooth function of temperature. Therefore, 
the numerical differentiation can be carried out without 
problem. 

The specific heat C and entropy S of ladders and 
tubes are plotted against T in Figs. 1 and 2, respec- 
tively. For both ladders and tubes, the overall behav- 
ior is not sensitive to q and is reminiscent of the two- 
dimensional ± J Ising model 7-10 except for the low tem- 
perature regime. For tubes, the entropy oscillates with 
the number of legs at low temperatures. This oscillation 
is not observed in ladders. This feature is also reflected 
in the low-temperature behavior of the specific heat as 



Fig. 8. Types of excitations: (a) droplet, (b) edge droplet, (c) 
domain wall, and (c') closed domain wall compatible with the pe- 
riodic boundary condition along the rungs. The spins in the shaded 
region are inverted relative to the ground-state configuration. 



shown in the inset of Fig. 1(b). 

To observe the g-dependence of the physical quantities 
in the low temperature limit clearly, the ground-state 
entropy is plotted against q in Fig. 3. The entropies of 
ladders and odd-leg tubes behave almost similarly, while 
the even-leg tubes have extra entropy in the ground state. 
The ground-state energy also shows a similar oscillation 
as shown in Figs. 4. These quantities arc plotted against 
\/q in Figs. 5 and 6. The results for ladders and tubes 
approach the values for the two-dimensional ±J Ising 
model 7 plotted by the open right-directed triangles with 
the increase of q. 

To observe the low-temperature asymptotic behavior 
of the specific heat, the quantity ln(CT 2 /7VJ 2 ) is plotted 
against J/T in Fig. 7(a) for ladders and (b) for tubes. 
The data for the ladders and the odd-leg tubes are close 
to each other and decrease with q. At low temperatures, 
they behave as CT 2 ~ exp(— 2 J/T) suggesting that the 
lowest excitation energy is 2J. The data for the even- 
leg tubes are below those for ladders and odd-leg tubes. 
This is consistent with the result that the ground-state 
entropies for the even-leg tubes are larger than those 
for other cases. For the even-leg tubes, the specific heat 
increases with q. At low temperatures, they behave as 
CT 2 ~ cxp(— 4J/T) suggesting that the lowest excita- 
tion energy is 4J. 

These features are understood by considering the ele- 
mentary excitations of the present models, which can be 
classified into the following three types: 

(a) droplet excitation, 

(b) edge droplet excitation, 

(c) domain wall excitation. 

Schematic pictures of these excitations are given in Fig. 
8. An excitation corresponds to the state with all spins in 
the shaded region inverted relative to the ground state. 
If an excitation has a vanishing excitation energy, it con- 
tributes to the ground-state degeneracy 

The excitation energy of a droplet excitation is a mul- 
tiple of 4J including zero, because its boundary always 
contains even number of bonds. The edge droplet exci- 
tation is allowed only for the ladders and its excitation 
energy is a multiple of 2J. The excitation energy of a 
domain wall is a multiple of 2J in the ladders. In the 
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tubes, however, only the closed domain wall is compati- 
ble with the periodic boundary condition along the rungs 
as shown in Fig. 8(c'). In this case, the excitation energy 
of a domain wall is a multiple of 4 J for the even- leg tubes 
and an odd-integer multiple of 2 J for the odd-leg tubes. 
Therefore, the lowest nonvanishing excitation energy is 
4J in even-leg tubes and 2 J otherwise. 

The above classification of excitations also helps to 
understand the excess ground-state entropy for the even- 
leg tubes. An even-leg tube can be formed by connecting 
two edges (a = 1 and q) of an even- leg ladder. If the even- 
leg ladder has a domain wall with energy 2 J, we can 
connect its both ends by inserting a vertical boundary 
between two edges. Thus, the domain wall in the even-leg 
ladder is converted into a closed domain wall in the even- 
leg tube whose excitation energies are multiples of AJ 
including zero. The contribution from these zero energy 
excitations can be interpreted as the excess entropy. It 
should be noted that this mechanism does not work for 
odd-leg cases, because the zero energy domain walls are 
not allowed in the odd-leg tubes. In the ladder, there is no 
constraint by the periodic boundary condition along the 
rungs. Hence, it is natural that the ground-state entropy 
and energy of the ladders vary smoothly with the number 
of legs. 

5. Summary and Discussion 

Finite temperature properties of the multileg ± J Ising 
ladders and tubes are investigated using the statistical 
transfer matrix method extending the method of Mat- 
tis and Paul. 4 It is found that the ground-state entropy 
shows an oscillating behavior with the number of legs in 
the tubes, while it decreases monotonically in the lad- 
ders. Corresponding behaviors of the specific heat and 
ground-state energy are found. From the numerical re- 
sults for the specific heat, it is found that the lowest ex- 
citation energy is 4 J for the even-leg tubes, while it is 2 J 
for other cases. The physical interpretation of these re- 
sults is given by analyzing the structure of excited states. 

The free energy of the two-leg ladder is calculated by 
MP using an approximate solution of their recursion rela- 
tion. In the course of the present investigation, however, 
we found that it substantially deviate from our numerical 
solution for ladders with q — 2. In addition, according to 
our numerical solution, the free energy does not approach 
the value of the ground-state energy predicted by MP. 
Actually, we found that the estimation of the ground- 
state energy by MP should be corrected. The corrected 
derivation of the ground-state energy and the numerical 
results for the temperature dependence of the free energy 
are given in Appendix. 

Recent investigations for the two-dimensional ± J Ising 
model suggests the power law temperature dependence 
of the specific heat in spite of the finite energy gap. 11 ' 12 
This anomalous behavior is attributed to the presence of 
the infinite rigid spin cluster with fractal dimension. In 
the finite width ladders and tubes, the power law behav- 
ior is excluded at low temperatures as shown in Fig. 7. 



However, this figure also shows that the specific heat 
crosses over from the high temperature regime, where the 
difference between the even-leg tubes and other cases is 
insignificant, to the low temperature regime, where this 
difference becomes significant. The crossover tempera- 
ture decreases with the increase of q suggesting the pos- 
sibility that it tends to zero in the limit of q — > oo. If this 
scenario is valid, the 'high temperature' regime can per- 
sist down to zero temperature in the limit q —¥ oo and the 
region with exponential temperature dependence would 
shrink to zero, allowing the power law behavior in the 
two-dimensional ± J Ising model. 

There are many possible extensions of the present 
model. In general, the magnitudes of the rung and leg 
interactions should be taken unequal. Similarly, the mag- 
nitudes of the ferromagnetic and antiferromagnetic inter- 
action should be unequal. The probability of each type 
of bonds can be different. The quantum effect would be 
most important in application to the real ladder and tube 
materials at low temperatures. The investigation of these 
effects on the present model is left for future studies. 
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Appendix: Ground State Energy of the ± J Ising 
Two-leg Ladder 

The ground-state energy of the ± J Ising ladder has 
been calculated by Derrida et al. 5 for q = 3 and Kad- 
owaki et al. 6 ' 13 for q — 2 and 3 using the zero tempera- 
ture transfer matrix method. Here, we present a simple 
derivation for the case q = 2 correcting the error of MP. 

When we assign ± J on the bonds of a ladder randomly, 
the probability that a plaquctte consisting of four spins 
on two neighbouring rungs is frustrated (or unfrustrated) 
is 1/2. If we denote the frustrated and unfrustrated pla- 
quette by F and U, respectively, each bond configura- 
tion is associated with a series of letters F and U, which 
can be identified by numbers of successive F and U as 
{^ui, n-pi, n\j2, riF2-"}- We assume the first plaquette is 
U without affecting the conclusion in the thermodynamic 
limit. We call a cluster of plaquettes consisting of succes- 
sive F's (U's) bounded by U (F) on both sides, a F(U)- 
cluster. 

In this representation, the ground-state energy for a 
bond configuration which corresponds to the sequence 
{n Vil riFi ■ i = 1, N c /2} is given by 



Nc/2 

E=J2 {Ev(nv,i) 



E F (n F ,i)} 



(A-l) 
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where A c is the number of clusters. Ep(n) and E\j(n) 
are the ground-state energies of F- and U-clusters with 
length n, respectively. For U-clusters, it is obvious that 



Eu(n) = -(3n + l)J. 



(A2) 



For F-clusters, it is energetically advantageous to put as 
many unsatisfied bonds as possible on the rungs which 
are shared by two F-clusters. This point was missed by 
MP. For even n, all unsatisfied bonds can be put on the 
rungs, while for odd n, one unsatisfied bond must be on 
a leg. Therefore, we find 

E F (n) = -{2n-l-mod(n,2)}J. (A-3) 

Note that the energies of the rungs on the boundaries 
between F- and U-clusters are counted in E\j. 
The total number of spins A is given by 

N c /2 

N = 2J2 (»ui + n Fi ) . (A-4) 

Since A is a macroscopic quantity, we can regard (A) as 
the actual total number of spins. The probability that a 
sequence of n letters appear is 1/2™. Hence, we find 



(nui) = (n Fi ) = (n) = ^ — = 2 



(A5) 



to obtain (A) = 4A C . Similarly, the ground-state energy 
per spin is calculated as, 

(E) = N C J 
A ~ 2A 
which is lower than the value —J predicted by MP. 



¥ = ^rH(») + M",!») 



7 4 (A-6) 
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To check the consistency with the finite temperature 
calculation, we have estimated the free energy for q = 2 
by iterating (15) and (16) for all possible realizations 
of {t n .a}- To make clear the comparison with Fig. 4 of 
MP, we plot the excess free energy AF = F — F pse udo 
for q = 2 against J/T in Fig. A-l, where i^pscudo is the 



pseudoanealed free energy 



F, 



pscudo 

A 



7 i -ln2 



In 2 cosh 



(A-7) 



defined by MP. The result substantially deviates from 
the approximate solution of MP. Also, as T tends to 0, 
it approaches J/3 as expected from eq. (A-6) instead of 
J/2 which is predicted by MP. 
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Fig. A-l. Excess free energy density AF/NJ (solid line) as a 
function of J/T calculated by direct numerical iteration of (15) 
and (16) for q = 2. The inset shows the low temperature behavior 
of AF/NJ plotted against T/J. The asymptotic approach to the 
exact ground-state value 1/3 (dotted line; right- directed triangle 
of inset) is observed. 



